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We investigate numerically the dynamics of two-dimensional Euler and ideal magnetohydrody- 
namics (MHD) flows in systems with a finite number of modes, up to 4096 2 , for which several 
quadratic invariants are preserved by the truncation and the statistical equilibria are known. Initial 
conditions are the Orszag-Tang vortex with a neutral X-point centered on a stagnation point of the 
velocity field in the large scales. In MHD, we observe that the total energy spectra at intermediate 
times and intermediate scales correspond to the interactions of eddies and waves, -&r(fc) ~ k~ 3 ^ 2 . 
Moreover, no dissipative range is visible neither for Euler nor for MHD in two dimensions; in the 
former case, this may be linked to the existence of a vanishing turbulent viscosity whereas in MHD, 
the numerical resolution employed may be insufficient. When imposing a uniform magnetic field to 
the flow, we observe a lack of saturation of the formation of small scales together with a significant 
slowing-down of their equilibration, with however a cut-off independent partial thermalization being 
reached at intermediate scales. 



PACS numbers: 52.65.KJ, 47.27.Ak,05.20. Jj 

I. INTRODUCTION 

A theory of turbulent flows is still eluding us, and yet 
such flows are ubiquitous in nature and are an integral 
part of the problem of weather prediction or of climate 
assessment, as well as for the formation and prediction 
of extreme events such as tornadoes and hurricanes; they 
are also an important element in the understanding of 
the dynamics of planets and the hcliosphere, stars and 
beyond. What is lacking is a statistical description of 
the small scales, and a prediction of long-time large-scale 
dynamics with ensuing modified transport properties. Of 
course, since Onsager (1949) pQ who considered an en- 
semble of point vortices in two space dimensions (2D), 
and T.D. Lee (1952) [5] who pioneered the study of the 
behavior of a truncated system of modes both for inviscid 
fluids and ideal magnetohydrodynamics (MHD) in three 
dimensions (3D), we know that equipartition of energy 
among the modes obtain in the simplest (3D, non heli- 
cal) case, whereas spectra peaking at the gravest mode 
in the 2D case are possible; these latter solutions can 
be viewed as precursors of inverse cascades toward large 
scales, observed in 2D for forced dissipative fluids and 
for MHD [3l |4j (see [5] for a review), as well as in 3D 
MHD j6] (see |7j for a review). Numerous extensions of 
the 2D problems introducing potential vorticity lead to 
the same prediction of an inverse cascade of energy, as in 
the quasi-geostrophic case. 

It has been argued that inviscid dynamics is in fact a 
good indicator of the behavior of turbulent flows in the 
presence of forcing and dissipation. This is certainly the 
case for scales larger than the forcing scale lp — l-njkp 
when an inverse cascade is present, but evidence for scales 
smaller than £p has been lacking until recently. It was 



shown, using high resolution ideal runs in 3D for Eu- 
ler flows that at intermediate times and intermediate 
scales, a Kolmogorov energy spectrum (hereafter, K41), 
E{k) ~ fc~ 5 / 3 , is observed, as well as a dissipative range 
at the end of the K41 spectrum; the latter range can be 
attributed to an eddy viscosity due to the thermalized 
high-k modes, E(k) ~ k? [S] (see also [5] in the heli- 
cal case when velocity-vorticity correlations are included, 
and [TU] in the 3D rotating case, with or without helicity) . 
One of the most striking result on the truncated Euler 
equations is thus that the thermalized zone progressively 
extends to large scales finally covering the whole spec- 
trum as expected, but that the intermediate scales are 
found to follow an effective Navier-Stoke equation [§1 E] 
even though dissipation is absent. 

Two-dimensional MHD is special for at least two rea- 
sons. First, the 2D ideal case possesses an infinite num- 
ber of invariants, even though only the quadratic ones 
are being preserved in the truncated system. Neverthe- 
less, the afore-mentioned predictions of inverse cascades 
have been verified in many instances [5]: one can indeed 
argue that other invariants, such as any power of the vor- 
ticity field integrated over space, will be dissipated effi- 
ciently in the realistic case in which forcing and viscosity 
are included. Furthermore, in MHD, there is a contro- 
versy as to what is the energy spectrum in the forced 
dissipative case: it could be a classical Kolmogorov spec- 
trum, either isotropic or anisotropic, in that latter case 
with a dependence on k± (referring to the direction per- 
pendicular to an imposed strong uniform magnetic field 
bo, of magnitude bo [H]). Or it could be an Iroshnikov- 
Kraichnan spectrum (hereafter, IK [13l [14]) stemming 
from the interactions of Alfven waves and turbulent ed- 
dies and leading to the slowing-down of the nonlinear 
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cascade when waves are strong; in fact, this latter so- 
lution is compatible with the prediction of weak turbu- 
lence theory for MHD [T5] when the resulting spectrum 
is isotropized. Moreover, it has been argued by several 
authors that there is no complete universality in MHD, 
including in the absence of forcing and with bo = (see 
[16j and references therein). 

In this paper, we thus investigate the properties of ideal 
two-dimensional MHD and Euler flows. The next sec- 
tion explains the procedure we follow and recalls several 
known properties of ideal 2D MHD; in we give re- 
sults in the absence of imposed magnetic field, bo = 0. 
In §IV[ we analyze the slowing-down of the dynamics for 
bo 7^ 0; finally, SjV] briefly examine the structures that 
develop in space, and |VI| is the conclusion. The special 
(b = 0) Euler case in two dimensions is treated in the 
Appendix. 



and for /(a) = 0,g(a) — —a; the former one, denoted 
A =< a 2 >, is the conserved square magnetic potential 
and the latter is called the cross helicity which can also 
be written as H c — J cPxh ■ u. As remarked in [5], the 
invariants in MHD do not go over smoothly into the con- 
served quantities of the hydrodynamics equations, except 
for total energy. We can expect then a different behavior, 
even if the magnetic field is weak initially. 

The truncated MHD equations for the pair of Fourier 
modes -0k and at, with k 6 [k m in,k max ], are defined in 
a similar way to the truncated Euler equation: 

dtp 1 

^^2(P X q^l^p^ - a P a q ]<Sk.p+q (7) 



dt 

dak 
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II. THE PROCEDURE 

We begin by writing the ideal MHD equations in the 
incompressible case (V ■ u — 0, where u is the velocity 
field) by introducing a pseudoscalar potential ip(x,y,t), 
the stream function, and the (scalar) magnetic potential 
a(x,y,t), with u = V x ip and b = V x a, b being the 
magnetic induction, also divergence-free (V • b = 0) in 
the absence of magnetic monopolcs: 



chp_ 
~di 
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di 



V 2 
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where {/, g} — d x fd y g — d x gd y f is the usual Poisson 
bracket. Regrouping some terms in equation ([2]), it can 
be shown that the magnetic potential also satisfies: 



da 

di 



u • Va = 



(3) 



and therefore a is advected as a passive scalar by the 
fluid, even though the Lorentz force j x b acts on the 
fluid and breaks the conservation of vorticity in 2D, (with 
j = V x b the current density) . Equations ([I] [2| conserve 
the total energy 



E = \ I |u|- 



(4) 
(5) 



they also have an infinite number of conserved quantities, 
the Casimirs, of the form 



C 



d 2 x[/(a) + VV<?(a)], 



(6) 



where / and g are arbitrary functions. Among them, two 
remarkable invariants are obtained for f(a) = a 2 ,g(a) = 



with 5k, r the Kronecker delta and with Fourier modes 
satisfying ipk — 0, at = if |k| > fc max ; for a computa- 
tional box of length 2tt, we have k m i n = 1, and with a 
de-aliasing using the usual 2/3 rule, k max = N/3 where 
N is the number of modes per dimension (we assume a 
box with a unit aspect ratio). 

This truncated system only conserves the quadratic 
invariants, which can be written in Fourier space as: 



E 
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(9) 
(10) 

(11) 



The absolute equilibrium is the equipartition distribu- 
tion of the invariants as first derived in [17] (see also [5]), 
with a, jH and 7 Lagrange multipliers, i.e. the parameters 
associated with the three invariants (note that He is not 
definite positive, and that 7 is a pseudo-scalar): 



aE + (3A + 7# c . 
This leads to the following equilibrium spectra: 

Eu(k) = ^[k 2 a + f3] , 



„ , , s 7rfc , 9 
Eb(k) = —k 2 a 



H c (k) 



2vrfc 



k 2 1 



(12) 

(13) 
(14) 
(15) 



with 



V = k 2 (a 2 - 7 2 ) + a/3 . 



For small values of j3, we have equipartition of kinetic 
and magnetic energy E u (k) = E^(k) ~ Trka/(a 2 — 7 2 ) at 
all scales. 
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III. THERMALIZATION IN THE ABSENCE OF 
A UNIFORM MAGNETIC FIELD 

A. The Orszag-Tang configuration in 2D MHD 

In order to study the thermalization of ideal MHD in 
two dimensions through a direct cascade of energy to 
small scales, we now resort to a numerical study. The 
code is a standard pseudo-spectral code with periodic 
boundary conditions; the temporal scheme is a Runge- 
Kutta time-stepping of fourth order that is known to ac- 
curately conserve the energy. The time steps used in our 
computations are given in Table|TJ note that they are sub- 
stantially lower than the CFL condition so that energy 
be preserved by the scheme for long times. 

We take as initial conditions the so-called Orszag-Tang 
(OT) vortex in 2D defined by the potentials [15] : 

t(j(x,y) — 2(cos k x + cos k y) , (16) 
a(x,y) = 2 cos k x + cos 2k y ; (17) 

with fco = 1 this vortex, made up of a neutral X-point 
for the magnetic field centered on the stagnation point of 
the velocity field, is concentrated in the large scales with 
E u = Et, = 2. at t = 0; hence, the initial rms values of 
the turbulent fields are u rms ~ 6 rms ~ 2. 

The short time dynamics is identical to the one stud- 
ied in jl9 at a resolution of N 2 = 512 2 Fourier modes 
and until the final time Tf — 1. Here, we integrate 
the truncated MHD equations until Tf — 20 using a 
grid of TV 2 = 2048 2 points. Note that the structures 
that develop at early times in physical space, namely a 
quadrupole in vorticity and a current dipole, were stud- 
ied in detail in [19 and are only briefly commented upon 
in^Vj 

The temporal evolution of the kinetic and magnetic 
energy is shown in Figjl] there is at first an energy ex- 
change between the kinetic and magnetic energy, with 
the magnetic energy dominating its kinetic counterpart. 
As thermalization is reached, after a time ~ 5, the ex- 
changes die out since the thermalized solution has no 
nonlinear energy flux; E b /E u settles at ~ 1.7, leading to 
Mrms ~ 1-5, &rms ~ 2.5. Such a moderate excess of mag- 
netic energy is often observed in the Solar Wind [20]; it 
corresponds to a slight departure from equipartition that 
can be attributed, in the context of ideal flows, as due to 
the effect of the conservation of the magnetic potential, 
see Eqs.( 13 14 1. We define as usual the total energy and 



total enstrophy spectra by summing the basic fields on 
circular shells of width Ak = 1: 

E T (k,t) = \ ]T [|u(k',t)| 2 + |b(k',i)| 2 ](18) 

fc_ <|k'|<fc + 

MM) = \ E [l^(k',t)| 2 + |j(k',i)| 2 ,(19) 

<|k'|<fc + 

with k± = k±Ak/2 ; u(k,t), w(k,t), b(k,t) and j(k,t) 
are the Fourier transforms of the velocity, vorticity, mag- 




FIG. 1: Top: Temporal evolution of kinetic and magnetic 
energy, E u (in blue) and Eb (in red) for a run at a resolution 
of 2048 2 grid points (Run 3, bo = 0); the dashed vertical line 
indicates the final time of the numerical integration reported 
in |19j on 512 2 grid points. Bottom: Temporal evolution of 
Et(k — 1)/E u (k — 1) (solid line) and of total enstrophies 
fit/fin =< j 2 > / < uj 2 > (dash line) for the same run; as 
long as the large-scale modes have not equilibrated, non-linear 
transfer to the magnetic energy persists whereas equipartition 
is reached rapidly in the small scales. 



netic field and of the current density, with j : z = V x he z . 
In Fig. [I] (bottom) , the ratio of the kinetic and magnetic 
energy spectra at k = 1 is displayed, as well as the ratio 
of < j 2 > / < oj 2 > (corresponding to the ratio of ki- 
netic and magnetic dissipation when equal viscosity and 
resistivity are re-introduced in the equations). Observe 
that magnetic energy dominates at large scale, a sign of 
a plausible inverse cascade associated with the magnetic 
potential and confirming that the defect in equipartition 
is indeed due to the large-scale behavior of the system, 
since almost exact equipartition is observed at small scale 
(see e.g. Figj3|. 

Concerning spectral variations, at early times as dis- 
played in FigJU the spectra show a clear kr 2 depen- 
dence, that can be associated with the formation of quasi- 
singular current and vorticity sheets |19) . This behavior 
is exhibited by all spectra (kinetic, magnetic and total 
energy). It is also present (data not shown) when a uni- 
form magnetic field is imposed to the flow, granted its 
magnitude is not too large compared to the rms values 
of the fluctuating field: for bo = 0.5, a similar spectrum 
obtains whereas for bo — 4, the spectrum is much steeper, 
due in part to a slower dynamical evolution (see next Sec- 
tion). This k~ 2 spectrum associated with sharp current 
and vorticity sheets persists as long as convergence of 
the partial differential equations is assured, i.e. until the 
thickness of current and vorticity structures is compara- 
ble (but still larger) than the grid size (see also the dis- 
cussion concerning the logarithmic decrement displayed 
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b = 2 


b = 4 


b = 8 
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b = 16 
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256 2 , Tf = 
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dt = 1/16000 
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512 2 , T f = 


60 


dt = 1/1600 


dt = 1/3200 


dt = 1/3200 


dt = 1/6400 


dt = 1/8000 


dt = 1/8000 
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2048 2 , Tf = 


= 20 


dt = 1/20000 
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2048 2 , Tf -- 


= 5 




dt = 1/10000 


dt = 1/20000 


dt = 1/40000 




dt = 1/40000 
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4096 2 , T f = 


1.2 


dt = 1/20000 













TABLE I: Time step dt for the various ideal MHD numerical simulations presented in this paper. In all cases, the total energy 
is conserved at better than 0.1% at the final time of the computation Tf. Note that b rms ~ 2.5 when bo = 0, so that all cases 
have a dominant imposed field bo, except the one with bo = 2 which can be viewed as an intermediate case. 



t=0.8 




FIG. 2: Kinetic (x, blue) and magnetic (circle, red) energy 
spectra (color online) at t ~ 0.83, in the early phase of the 
flow for Run 3 with bo = 0. The black dashed line corresponds 
to fc~ 2 , a spectrum expected in the presence of quasi-singular 
structures. Note that magnetic energy dominates at all scales 
at that time. 

in Fig. |). 

It was shown in [8] in the context of three-dimensional 
(D=3) Euler flows that, as time evolves, another dy- 
namics takes place: once the flow behaves as a trun- 
cated system of modes and thermalization begins, the 
modes following a fc D_1 law at small scale for the three- 
dimensional Euler equations produce an eddy viscosity 
for larger scales; a quasi-turbulent regime follows at in- 
termediate scales and intermediate times (i.e. before 
thermalization occurs everywhere), with an inertial range 
close to a Kolmogorov law corresponding to the forced 
dissipative case. 

Similarly in the present 2D case, once the smallest re- 
solved scale has been reached, the overall solution be- 
comes noisier, the vorticity and current sheets curve and 
interact, and the Fourier spectra evolve as well (struc- 
tures are displayed in Fig{9] below) . 

The kinetic and magnetic energy spectra, compensated 
by a fc 3 / 2 power law, are shown in Figj3]for several times, 
in the range t ~ 2.5 to t ~ 20. As in the case of ther- 
malization for the truncated Euler equation evolution in 
hydrodynamics with large-scale initial conditions, a clear 
scale separation now appears in these spectra. 

It is possible then to define a wavenumber ktt(t) where 



the thermalized fc 1 power-law zone starts (numerically 
determined by seeking the minimum of E{k)). The total 
thermalized energy and enstrophy are thus defined as: 

Ethif) = £ E T (k,t), Q th (t) = ]T fir(M). (20) 

We now can estimate the effective energy dissipation rate 
as e th (t) = Similar to the 3D Euler case [El M, this 

quantity presents a maximum, here for t f=a 2.5 (data not 
shown). Note that this time also corresponds to the in- 
terval of time during which the total dissipation, when 
viscosity and magnetic resistivity are added to the primi- 
tive equations, is quasi-stationary (see [3T]). The scaling 
laws for a K41 spectrum (dotted line), and for the IK 
spectrum (dashed line) defined respectively as: 

E% 41 (fc) « E™ (fc) ~ C K e 2/3 fc- 5/3 , (21) 

and 

Ei K (k) a El K {k) = C IK [sb lms ] 1/2 k~ 3/2 (22) 

are also displayed in Fig. [3]at large scales, with 6 rms ~ 2.5 
the root mean square value of the fluctuating magnetic 
field. The numerical constants Ck and Cjk appearing 
in front of the K41 and IK spectra are evaluated (see 
Fig. [3ja) by fitting the large-scale part of the spectra at 
an early time after the first Alfvenic energetic exchange, 
t ~ 2.5 (note that the turbulence spectrum is already 
visible at t ~ 2.0, not shown). We find for these constants 
Ck = 2 and Cik = 0.8 at the time of the maximum of 
the effective dissipation £th(t). 

Therefore, we can safely conclude, on the basis of the 
examination of the results plotted in Fig. [3] that the IK 
spectrum is obtained in the intermediate spatio-temporal 
range of ideal dynamics in MHD in two dimensions for the 
OT vortex. Also note the good equipartition of kinetic 
and magnetic energy in the thermalized range, and the 
domination of magnetic energy at large scale, again a 
sign of a plausible inverse cascade associated with the 
magnetic potential. 

B. Random initial conditions 

Similar results obtain for random initial conditions, as 
can be seen in Fig(4] plotting spectra in that case for 
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FIG. 3: Energy spectra for various times, compensated by fc 3//2 : kinetic (x, blue) and magnetic (circle, red). Dashed and dotted 
lines represent fit for the K41 and IK solutions respectively, the constants appearing in these spectra being fixed at the earliest 
time shown, t ~ 2.3. Grid of 2048 2 points, for Run 3 with bo = (see Table III. Note the fast equipartition at small scale, the 
IK spectrum for more than one decade at intermediate scales and intermediate times, and the persistence for long times of the 
domination of magnetic energy at large scale. 



a normalized global velocity-magnetic field correlation 
(H c /E) that takes values of (left) and 0.8 (right): a 
thermalized spectrum is observed at small scale, and a 
turbulent IK-spectrum at large-scale. However, note that 
non-universality has been obtained in MHD |16| in the 
decaying dissipative K41 spectrum was also pre- 

viously obtained in decaying 2D MHD turbulence with 
random initial conditions using a Lagrangian model |22| . 



C. The pseudo-dissipative range, or lack thereof 

Remark that the end of the pseudo-inertial range is not 
followed by a sharp decrease in energy before the ther- 
malized Er{k) ~ k small-scale spectrum, contrary to the 
3D Euler case [9]. However, computing a dissipation 
wavenumber using the IK spectrum and the correspond- 
ing evaluation of the transfer time that leads to the IK 
spectrum, gives idiss ~ v 2 ^ z , where in our ideal case the 
viscosity is replaced by a turbulent expression based on 
the thermalized energy, v turb ~ \fE[h/k max ~ l/k max . 
Similarly, one can recall that the eddy viscosity computed 
with the EDQNM closure gives a non-zero contribution 
23J. However, several remarks are in order. First of all, 
the small-scale velocity leads to no contribution to an 
eddy viscosity for the dynamics of the large-scale velocity 



field [5] , so the sole contribution to dissipation of E u will 
stem from the small-scale magnetic field. Furthermore, 
the eddy resistivity contributions of the small-scale ve- 
locity and magnetic field exactly compensate each other 
when E u = Ef, in the small scales, which is the case here 
(with small /?). However, one can compute the correc- 
tion to equipartition which is known to follow a k~ 2 law 
in the dissipative case [23]; therefore, one could expect a 
non-zero contribution to turbulent viscosity in 2D MHD 
as well. Finally, another argument can be put forward 
to explain the lack of sharp decrease of the spectra be- 
fore the thermalized range, namely that there are not 
enough modes in quasi-equilibrium to produce a suffi- 
cient amount of effective dissipation. Indeed, in the 3D 
case, it was shown that a number of roughly 256 3 modes 
was necessary to see this internal decrease of the energy 
spectrum before the thermalized spectrum; this would 
correspond to a computation on a grid of 4092 2 points 
in 2D. However, observe that using a resolution of 2048 2 
points, a slight dissipative zone seems to appear at t = 2.3 
(for k fa 50). This point will await further study. 

To try to understand further the lack of dissipation 
range, we also resorted to an examination of the two- 
dimensional Euler case that is given in some detail in 
Appendix [X] The truncated Euler equations relax to- 
ward the statistical equilibrium in an analogous way to 
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FIG. 4: Energy spectra for random initial conditions with normalized velocity-magnetic field correlation of (left) and 0.8 
(right), 3/2-compensated kinetic (x, blue) and magnetic (circle, red). The dotted lines correspond to an IK spectrum, with 
Cik = 0.8 (left) and Cm = 0.5 (right), see eq. (221. 



the three dimensional case. The main difference is the 
presence of a direct cascade of enstrophy. This quan- 
tity plays the role of the energy in 3D, thermalizing in 
equipartition at large wave- number and yielding a k~ l 
law in the inertial zone. Again, a remarkable difference 
with the 3D Euler case is the absence of a dissipative 
zone that is due here to a vanishing 2D eddy viscosity. 



IV. DYNAMICAL SLOWING DOWN IN THE 
PRESENCE OF A UNIFORM MAGNETIC FIELD 

In the presence of a strong imposed uniform magnetic 
field of amplitude 6q , it is known that the dynamics is 
slowed down, including in the ideal case |19) ; this is in fact 
at the basis of the argument of Iroshnikov and Kraich- 
nan for an energy spectrum in MHD different from the 
Kolmogorov spectrum, and it is the feature on which the 
weak turbulence development for MHD relies upon [T5] . 
using the smallness of the ratio of the Alfven time to the 
eddy turnover time 



Rwt = ta/tnl , 



(23) 



with ta = L /bo and t^l = L /Uo, L and Uq being 
respectively the characteristic large scale and velocity of 
the flow. It is also claimed in [T5], in the framework of 
ideal 2D MHD, that in fact the development of small 
scales in the presence of a sufficiently strong b is ar- 
rested, with a smallest excited scale £b , that depends on 
bo and that can be larger than the smallest resolved scale 
of the flow in a computation at a given resolution. 

In view of the increased power of computers available 
today, we revisit the effect a uniform held has on the 
formation of small scales in the ideal case of 2D MHD. 
Compared to the work in |19j . we are now performing 
computations for longer times, for different values of the 
magnetic field bo and for higher resolutions, using here 
grids up to N 2 — 4096 2 points. Note that, for a magnetic 
field aligned with the s-axis, eqs.([Tj2| are modified by 
performing the substitution a — > a + b y. Remark that 
now the square magnetic potential A =< a 2 > is no 



longer conserved and the absolute equilibrium spectra 
given m Eqs.( pp5| ) are obtained in this case by setting 
= 0. 

In Fig(5]we give the temporal evolution of two modes, 
one in the middle of the resolved range (k = 13), and the 
other one at the end [k = k max — N/3), both normalized 
by E(k=l,t=0), and for several values of the imposed field 
(see caption). We observe a delay in the early dynamics 
of the modes as bo increases, followed at long times by 
a saturation once equilibrium is reached (particularly so 
for E(k max )). 

The logarithmic decrement technique |25j is now ap- 
plied to quantify further the delay of the onset of the 
evolution as the amplitude of the imposed uniform field 
is increased. If the fields are regular, then the energy 
spectra must decay at least exponentially at large wave- 
number k. Based on this assumption, the logarithmic 
decrement S(t) is defined by the large k asymptotic of 
the energy spectra: 



E T (k) = c{t)k~ m{t) e~ 2mk 



(24) 



S(t) is measured by fitting the long wave-number range 
and the minimum admissible value is determined by the 
relation S(t)k max — 2 (S =twice the mesh). The tempo- 
ral evolution of S(t) for different values of the imposed 
magnetic field is displayed in Fig|6]in a log-lin plot. 

The presence of a strong magnetic field slows down the 
nonlinear interactions and S(t) remains above the min- 
imum admissible value for a longer time the stronger 
the value of bo, but it finally reaches the condition 
S{t)k max = 2 near t = 3 for bo = 4 (and t ~ 5.8 for 
b = 8, see Fig(6]). From this study, we can deduce 
that the system can reach statistical equilibrium even in 
the presence of a strong imposed magnetic field, but that 
the convergence toward such a state is considerably ham- 
pered. In Fig{6] (bottom), we also display the temporal 
evolution of the kurtosis for the current and vorticity and 
for several imposed uniform magnetic fields. Again, the 
delay in the formation of small scales is observed, as well 
as a tendency towards Gaussianity once the thermaliza- 
tion takes place. 
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FIG. 5: Temporal evolution of total (kinetic plus magnetic) 
energy in mode k = 13 (top) and k = k max (bottom) for sev- 
eral imposed mean fields: bo — (x), bo = 2 (circle), bo = 4 
(*), bo = 8 (triangle) and bo = 16 (star); both modes are nor- 
malized by the initial value of the energy in the gravest mode. 
Note the substantial delays in the evolution as bo increases, 
and the final saturation level. Grids of 512 2 points (Runs 2). 



It may appear somewhat surprising that a different 
conclusion is reached with the present data from what is 
argued in |19j on the basis of a quasi-regular behavior of 
MHD in two-dimensions in the presence of a large-scale 
magnetic field: indeed, such a field retards the nonlinear 
dynamics and can be seen as a bath of weakly interacting 
Alfven waves with a spectrum that can be derived ana- 
lytically in the case of a strong enough b , using weak 
turbulence theory |15) . However, it is well-known that 
the weak turbulence approach is non-uniform in scale: 
the small parameter of the problem, Rwt, should be 
evaluated taking into account that the eddy turn-over 
time gets smaller as smaller scales are reached, whereas 
the length scale of the imposed field remains infinite, by 
construction. Working out this relation in the case of 
the IK spectrum, E(k) ~ (e^o) 1 ^ 2 ^ -3 ^ 2 leads to a scale 
£b ~ 1 /bo with Rwt(^b) — 1, scale beyond which a 
classical small scale turbulent spectrum will develop; it 
should be noted, however, that the resolution of the com- 
putation must be such that £g is reachable accurately 
(2^sfc max > 1) in order to observe this phenomenon. 

We now examine numerically the spectral relaxation 
to equilibrium; several times are displayed in FigjTJ with 
dashed lines corresponding to a k +1 scaling. We first 
observe that, at early times, the spectra coincide at large 
scales, but that small-scale thermalization is delayed for 
stronger imposed fields. 



FIG. 6: Top: Temporal evolution of the logarithmic decre- 
ment 5(t) with bo = (x), bo = 2 (circle), bo = 4 (star), 
bo = 8 (triangle) and bo = 16 (*). Resolution of 2048 2 grid 
points, Runs 3 and 4. The horizontal solid line gives the small- 
est admissible value of <5 for the PDEs to be well resolved for 
these grids, and the vertical dashed line represents the time at 
which integration was stopped in [19] for a computation on a 
grid of 256 2 points. No saturation of the temporal evolution 
of the logarithmic decrement 8(t) is observed provided the 
numerical resolution allows one to compute for long enough 
times. Bottom: Kurtosis 72 = K.4/K2 (with Ki the cumulant 
of order i) for the same runs as above, of the current j and 
vorticity uj for different values of the background constant 
magnetic field; the current is shown with a solid line (and 
darker color), whereas the vorticity is given with a dash line 
(and lighter color, see left inset). In the right inset is shown 
the same plot in log — lin coordinates; 



At intermediate times, there is a domain of wavenum- 
bers in which the dynamics differs quite substantially ac- 
cording to the value of b . It appears that there is now 
a relaxation to equilibrium at intermediate wavenumbers 
(the k 1 scaling law seems to be followed, see the figure) 
before a sharp plunge in the spectrum due to the fact that 
small scales are not reached yet because of the substantial 
slowing down of the dynamics. This is particularly strik- 
ing at t ~ 5 for which the run with the largest value of 
the imposed mean field already presents a partial ther- 
malization at intermediate wavenumbers, as if the dy- 
namics was indeed seeing a truncation at wavenumbers 
much lower than the actual k max of the run, whereas for 
bo = 4, energy continues to flow to smaller scales. At a 
later time (t ~ 15), the bo = 4 case now may be seeing 
a pseudo-truncation in scale and partially thermalizes at 
intermediate scales, but it still cannot reach the smallest 
scales numerically available to the run. At the final time 




FIG. 7: Total energy spectra with bo = (x), bo = 2 (circle), bo = 4 (*), bo = 8 (triangle) and bo — 16 (star) at t = 1.2 and 
i = 5 (top, with a grid of 2048 2 points, Runs 3 and 4), and t = 15 and t = 60 (bottom, grid of 512 2 points, Runs 2). The 
dashed lines correspond to a k 1 scaling. Thermalization may occur for very long times, but the spectral index at intermediate 
scales/times for higher bo is far from the predicted value, with possibly a partial thermalization at these intermediate scales 
(see top, right). 



of the computation, t ~ 60, the three runs start to dis- 
play thermalization in a broader range of scales, although 
the evolution of the gravest mode is quite retarded again 
compared to the weak field case. However note that the 
exponent of the thermalized zone clearly depends on bo, 
(see Figj8ja); for b = 16, the exponent is close to 
indicative of a partial one-dimensional thermalization. 

FigurejHJb displays the time T* at which the small- 
est effective excited scale in the flow reaches the grid size 
(S(T t ) = 2/k max ). Observe that T* is an increasing func- 
tion of bo and no sign of saturation of this slowing-down 
is observed yet, although one might want to test higher 
values of bo as well. Including all points, a quadratic vari- 
ation of T* with bo is plausible. Note however, that the 
close to zero value of the spectral index for bo = 16 (sec 
Figj8]a) is in apparently contradiction with the expected 
k +1 power-law of equipartition of energy in 2D. It re- 
quires substantially longer times of integration to check 
whether the convergence toward the predicted statistical 
equilibria is particularly slow or whether another solution 
is obtained (as in the Fermi-Pasta-Ulam-Tsingou prob- 
lem [55]). In order to investigate this point, we performed 
an integration until t = 475 at moderate resolution. The 
spectral index for this run is displayed in Fig|8]c and it 
is found to asymptotically approach to 1, the predicted 
value of the 2D absolute equilibrium. 



V. STRUCTURES 

The structures that develop in the flow we study here 
have been examined in detail in |19| at early time. Here, 
we pursue this study at higher resolution and examine 
how these structures change with bo and with time. The 
structures are shown in Fig. [9] using grids of 2048 2 and 
4096 2 points (note that the color bars are described in the 
caption). At early times, the flows are made up of quasi- 
singular structures which persist until t < 1 for bo = 
(Figs. [9] top) and for longer times when the magnitude 
of the external field is increased (see Figs. [9] bottom, 
obtained for the case 60 = 4). Once the thickness of 
these structures reach the grid size, noise steps in but 
the structures continue to evolve until the system is able 
to reach complete thermalization. 

In the presence of an imposed field, when strong 
enough, structures align themselves in its direction. Fur- 
thermore, the appearance of noise is clearly delayed; at 
a given time, it is not yet visible as the amplitude of bo 
increases (contrast the current at t = 1.1 for 60 = (top- 
right) and bo = 4 (bottom-left)). Finally, note that in 
all cases with strong imposed field, the small-scale struc- 
tures display an intense folding and piling-up of sheets of 
opposite signs. 
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FIG. 8: Variations with imposed mean field bo- a) Spectral 
index a = — m (see Eq.(|24|) with S(t) = 0) of the thermalized 



-m (see Eq.p4| with S(t) 
part (k > 20) of the total energy spectra at t = 60 (see Fig[7|) 
for several bo- Orszag-Tang vortex, resolution of 512 2 grid 
points, Runs 2 (see Table b) Time T* as a function of bo; 
T* is defined as the time at which the grid size is reached by 
the computation, as measured by the logarithmic decrement 
S. c) Temporal evolution of spectral index for k > 20 with an 
imposed uniform field bo = 16 at resolution 256 2 (Run 1); note 
the slow evolution toward a = 1 as expected at equilibrium. 



VI. CONCLUSIONS 

We have investigated in this paper the dynamics of 
an ideal two-dimensional fluid in the MHD limit, in the 
presence or not of a uniform magnetic field bo, and 
we have shown the link to the dissipative driven (DD) 
case. For bo = and at intermediate times and in- 
termediate scales, a behavior observed in the DD case 
obtains, namely that the energy spectrum is that pro- 
posed by Iroshnikov and Kraichnan, with a /c~ 3 / 2 power 
law [T31 H3] for the Orszag-Tang vortex and for random 
flows, with in both cases initial conditions centered in 
the large scales. Furthermore, as already found in [19] 
but for shorter evolution times and lower resolutions, the 
formation of small scales is inhibited when bo =/= 0. 

Note that small-scale initial data was also studied in 
[27j where it was shown that the system reaches equi- 
librium by an eddy-noise mechanism. No inverse cascade 
was observed as mechanism of thermalization. In the case 
of MHD an equipartition of magnetic potential leading 
to a k 3 scaling-law for the magnetic energy spectrum and 
k 1 equipartition for the kinetic energy spectrum was ob- 
tained at large times. 

Extension of this work to the three-dimensional case 
in MHD may be of use for at least three reasons: 
(i) It has been shown in [TB] that, for initial conditions 






FIG. 9: The color bars (not shown) are such that red and 
magenta represent extreme values, green and blue moderate 
to low values of the fields: stream function and vorticity vary 
from -4 (red) to 4 (purple), with in the former case most of 
the points at values close to zero, whereas the vorticity in 
this range has all values with close to equal probability; the 
magnetic potential varies from -3, red, to 3, purple, and the 
current from -6, red, to 6, purple. Top: Current for bo = 0, 
t = 0.9 (left) and t ~ 1.1 (right). Resolution of 4096 2 grid 
points. Bottom: Current for bo = 4, t — 1.1 (left) and t ~ 3.7 
(right). Resolution of 2048 2 grid points. In the presence of 
even a moderate uniform field (bo/b rms ~ 1.6 here), the onset 
of thermalization that destroys structures is delayed and the 
flows align themselves with bo. 



that are identical from the point of view of the statistics 
(same energy, same velocity field, same equipartition be- 
tween kinetic and magnetic energy at t — 0, same initial 
conditions centered in the large-scales, same total mag- 
netic helicity, and with total velocity-magnetic field cor- 
relation between and 4%, in normalized value), three 
different energy spectra could emerge in the absence of 
imposed uniform magnetic field and forcing (decay case 
with non-zero viscosity and unit magnetic Prandtl num- 
ber) when considering three different initial conditions for 
the induction. Would they be observed as well at inter- 
mediate times- intermediate scales in the ideal case? It is 
plausible to think so, since a Kolmogorov spectrum is ob- 
served in the ideal 3D case for neutral fluids, and we ob- 
serve the Iroshnikov-Kraichnan spectrum in the present 
work, but it would be of interest to verify the lack of 
universality in MHD in the ideal case as well. 

(ii) Would a fast-decreasing spectrum at the end of the 
inertial range and before the thermalized range obtain 
in three-dimensional MHD, as it does in the 3D neutral 
case? It is argued in this paper that the lack of such a 
range is probably due to an insufficient number of ther- 
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FIG. 10: Temporal evolution of enstrophy spectra k 2 E(k) for 
2D Euler, at t = 14 (top) and t = 40 (bottom). Solid and 
dashed lines respectively indicate the enstrophy cascade (& _1 - 
scaling) and thermalized enstrophy at small scales. Again 
note, as in the 2D MHD case, the absence of a pseudo- 
dissipative range before equipartition sets in. 



tional Science Foundation. 

Appendix A: The fluid case in 2D 

We take here as initial conditions, b = 0, and for the 
stream function we have: 



ip(x, y) = — sin k a x sin k a y + — cos k^x, (Al) 

K a Kb 

with the parameters set to k a =l and kb~2. The tempo- 
ral evolution of the enstrophy spectra is displayed in Fig. 
[TO] As in 3D, a clear scale separation also appears: there 
is a progressive thermalization starting from the smallest 
scales, with the energy (vs. enstrophy) cascading to the 
larger (vs. smaller) scales. Note that the thermalized en- 
strophy increases from zero at early times to an amount 
of the order of the total enstrophy available in the sys- 
tem. Defining k t h as previously, and using the values of 
Qth and E^, we can compute the parameters (Lagrange 
multipliers) a and /3 from: 



E(k) = —— , fi(fe) = kE(k) . (A2) 



malized modes because the total number of modes in 
two dimensions is not very large, compared to the three- 
dimensional case at the same linear resolution; thus, an 
effective eddy viscosity does not obtain here, and such a 
3D computation would provide a test of this idea, 
(iii) Finally, the case of magnetic helicity Hm =< arj >, 
an invariant in ideal MHD in three dimensions, deserves 
a separate study. Indeed, it was shown in |24j that, per- 
forming a standard phenomenology a la Kolmogorov on 
Hm, the spectrum becomes HM(k) ~ fc -7 / 3 but recent 
studies have shown that different spectra obtain [28H30] 
with HmQz) ~ k~ 3 or steeper. The origin of this discrep- 
ancy is not completely understood; it could be related 
to a detailed equipartition between kinetic and magnetic 
modes in the energy and helicity parts of their spectral 
correlation functions. Again, an ideal study may help un- 
ravel the mechanisms at play in the dynamical evolution 
of MHD turbulence. 
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These Kraichnan absolute equilibria are displayed as 
solid lines at small scale in Fig. 10 (notice that they 



curve down at the larger scale of the thermalized zone for 
enstrophy); these solutions correspond to the /^-scaling 
of the high enstrophy containing absolute equilibria. 

The good agreement shows that the evaluation of k t h 
and the energy and enstrophy at the wavenumber de- 
scribe reasonably correctly the temporal behavior of the 
flow (in the latter case, as r^l^k -1 ). The enstrophy dissi- 
pation rate 7y t h can be estimated, defining it as the time 
derivative of fi t h- The respective spectra are also dis- 
played as solid lines at large scales in Fig. 
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that in the inertial zone both scaling law and prefactor 
are in good agreement with f2(fc) ~ 77?/ 3 A; — 1 . 

The relaxation dynamics of two-dimensional Euler 
turbulence is, mutatis mutandis, similar to the three- 
dimensional Euler case: a direct cascade of enstrophy 
(energy in 3D) followed by an equilibration of enstrophy 
(energy in 3D) at small scale. It can be shown [57] that 
the dissipation wavenumber estimated from the equilib- 
ria is of the order of the maximum wavenumber, thereby 
explaining in this case the absence of dissipation range 
as a buffer zone between the inertial range and the ther- 
malized range, contrary to the 2D MHD case. 
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